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We study two-body correlations in systems of identical bosons. We use a Faddeev type of de- 
composition of the wave function where all pairs of particles are treated equally. We focus on a 
new multi-particle Efimov effect at large scattering length, where infinitely many loosely bound 
many-body states appear. A confining external trap only allows a finite number of such spatially 
extended negative energy states inside the trap. The stability of a Bose condensate is determined 
by the decay into these model independent intermediate states which in turn decay into dimers. 
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Introduction. The novel theoretical formulation in 
was constructed to describe correlations in boson systems 
and applied for a realistic short-range repulsive interac- 
tion to Bose-Einstein condensates. The method goes be- 
yond the mean-field approximation J^, |j| and is as well 
applicable to attractive finite range potentials with very 
large scattering lengths where the Efimov effect occurs 

|^| . Experimental properties are available for various 
condensates from the first observations of both effectively 
repulsive || and attractive |7[| interatomic interactions to 
the recent observations jq, [Jof condensates in a magnetic 
field used to tune the effective interaction via a Feshbach 
resonance to almost any scattering length. 

The condensate § being an excited state of the full 
many-body system is clearly unstable. The three-body 
recombination into bound dimers is a dominating decay 
channel Jl(], O, [l2|. In the condensate this occurs in- 
dependently for neighbouring pairs and much more fre- 
quently in a coherent process best described as a collec- 
tive or macroscopic collapse [hi]. Two-body corre- 
lations therefore must be crucial for this collapse, which 
becomes more likely and eventually inevitable as the scat- 
tering length is increased. Unfortunately a theoretical 
description is hindered by the difficulties, especially pro- 
nounced for large scattering lengths, of finding the deci- 
sive correlated structure |2| [l5[ [l6|. 

A promising form of a correlated wave function sug- 
gested for nucleons |17 was recently extended to more 
general systems [Q. Another related formulation uses 
generalized hypcrspherical coordinates and an adiabatic 
expansion with the hyperradius as the adiabatic coordi- 
nate ]l9| . It was applied for many identical bosons, still 
only with a zero-range interaction and the lowest (con- 
stant and thus non-correlated) hyperspherical angular 
wave function. These crude approximations are removed 
in a novel method still using hyperspherical coordinates 
and adiabatic expansion, but now the wave function is 
determined from a variationally established equation [Q . 
The purpose of this letter is to investigate the structure 
of boson systems as a function of (large) scattering length 
for attractive finite range potentials with emphasis on an 
emerging novel many-body Efimov effect. 



Theory. A system of N identical particles of masses 
m, trapped in an external field approximated by a har- 
monic oscillator potential of angular frequency w, can in 
the center of mass frame be described by hyperspheri- 
cal coordinates, i.e. 3A^ — 4 hyperangles |TJ|, ^| denoted 
collectively by f2 and one hyperradius, p, given by (!]] 
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where fj are single-particle coordinates, R is the center 
< tJ — \fj — fi\ = \/2p sin cvjj with 
an varying between and 7r/2. The hamiltonian is 



of mass coordinate, and r. 
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where V is the two-body interaction. It separates into a 
center of mass part, a radial part, and an angular part 
depending respectively on R, p, and £1 m: 
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with H c 



= pj i /(2Nm) + Nmuj 2 R 2 /2, H p = T p + 
and h 2 h n /(2mp 2 ) = f a + J2i<j Vij, wh ere f p 
and Tq are related kinetic energy operators. A suitable 
expansion of the wave function is 
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where is an eigenfunction of the angular part of the 
hamiltonian with an eigenvalue h 2 X n (p) / (2mp 2 ): 



hfi^ n { Pi n) = \ n (p)$ n (p, O) 
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Neglecting couplings between the different n-channels 
yields the radial eigenvalue equation for the energy E n : 
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The second term in the radial potential U is a generalized 
centrifugal barrier. We now decompose the angular wave 
function $ in the symmetric Faddeev components cj> 
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where the restricting assumption is that the interparticle 
potentials only act in s-waves leaving only the depen- 
dence on the distance r.y = \/2p sin otij . The capability 
of this decomposition for large scattering length has been 
demonstrated for N = 3 by the proper description of the 
intricate Efimov effect [|[ 

The eigenvalue equation eq. (H) can by a varia- 
tional technique be rewritten as a second order integro- 
differential equation in the variable a 12 For atomic 
condensates the interaction range is very short compared 
to the spatial extension of the TV-body system. Then this 
equation simplifies even further to contain at most one- 
dimensional integrals. The principal interaction depen- 
dence is through the parameter as = m J °° drr 2 V(r)/h 2 , 
which is the Born-approximation to the s-wave scattering 
length a s . The validity of our approximations only relies 
on a small range of the potential whereas the scattering 
length can be as large as desired. 

Angular potentials. The crucial ingredient is the an- 
gular eigenvalue A obtained from eq. (0). We choose gaus- 
sian two-body potentials V(r) = Vo exp(— r 2 /b 2 ) of range 
b and strength Vo. We measure lengths in units of b and 
energies in units of h 2 / (2mb 2 ). We solve the equation by 
choosing discrete sets of mesh points and thus construct- 
ing a matrix equation ready for diagonalization Q . This 
computational scheme increases strongly with N and we 
therefore illustrate by relatively small values of N. 

Fig. [I] shows the angular eigenvalues for different 
strengths expressed in terms of a B and a s . For short 
range potentials finite at the origin all A approach their 
values at p — as p 2 . This is a model dependent region 
and therefore not considered in the following. When the 
potential is too weak to support a two-body bound state, 
(a B /b > —1.18934), the lowest Ao is negative and ap- 
proaches zero as a s / p with increasing p. In fig. |l| we also 
show the eigenvalue Ai corresponding to the first excited 
angular state. It is positive and approaches the value 
4(4 + 3iV - 5) as -a s /p for large p @ §l). 

Increasing the attraction of the potential a two-body 
bound state of energy B2 appears (positive a s ) and the 
corresponding lowest angular eigenvalue diverges propor- 
tional to 2mB2p 2 /h 2 when p — ► 00. The many-body sys- 
tem would obviously then prefer to recombine into bound 
dimer states. The first excited angular eigenvalue, sim- 
ilar to the dotted curve in fig. [I], is then responsible for 
the unstable state referred to as the condensate. 

At the two-body threshold the scattering length is in- 
finite. Then the lowest A approaches a negative constant 
Aoc for large p (p > pi, = 45^/Nb), see fig. |l|. Using 



N = 3, 10, 20 we find numerically Aoo(iV) w 24(A - 2) - 
5N(N — 1), where the form reflects our expectations of 
iV-dependence of kinetic and potential energies. 



-1000 



-2000 



-3000 



Ao, (-1.11,-15.7) 
Ax, (-1.11,-15.7) 
Ao, (-1.1893, 00) 
A , (-1.22,44.5) 
1 I 1 1 1 I 1 



200 



400 



600 



800 



1000 



p/b 



FIG. 1: Angular eigenvalues as functions of p for N — 20 
and various potential strengths given by (a B /b,a s /b). The 
long-dashed and the dotted curves correspond to the two low- 
est eigenvalues when the two-body potential is too weak to 
support a bound state. The potential giving rise to the short- 
dashed curve has only one bound two-body state and the solid 
curve corresponds to the two-body threshold of zero energy. 

These angular eigenvalues behave qualitatively differ- 
ent from those obtained in [|l9), where the expectation 
value of a (5-interaction of strength proportional to a s is 
computed with a constant angular wave function without 
any correlations. The eigenvalues are then necessarily 
proportional to a s / p for all p, i.e. always diverging when 
p — ► and converging to zero for p — ► 00. The differences 
are especially pronounced for large scattering lengths and 
when a bound two-body state is present. However, even 
in the qualitatively similar case in fig. |of a s /b= -15.7 
our A is lower than that of by a factor between 1.5 
and 2 for p/b G [10 3 ,10 5 ]. 

At large p, <& approaches a non-correlated angular wave 
function. Using a constant as in |19 the expectation 
value of our gaussian interaction gives for large p and N 
that A = A n _ c (p) = kN 7 / 2 a B /p with k = (3/2) 3 / 2 /V2^, 
i.e. we get ae instead of a s as assumed in For large 
a s our computed A bends upwards from the plateau of 
Aoo around a point p = pth determined by the large but 
finite value of a s . If p > pth then A can be estimated as 
An-c(p) with a s substituted for ob- The transition point 
is then given by Aoo = A n _ c (pth) which for large N yields 
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where the more realistic estimate for attractive potentials 
is indicated by exchanging a B — > a s , sec |l^, |2l| . 

Radial potentials. In eq. (Q) the external harmonic 
trap of angular frequency u> corresponds to a trap length 
b t = ^/ h/ (mu>) . We model the experimentally studied 
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systems H |9| of 85 Rb and 87 Rb-atoms with oscillator 
frequencies v — uj /(2tt) = 205 Hz and 200 Hz and inter- 
action range b — 10 a.u., thus yielding b t /b = 1442. 

We show in fig. |^ the radial potentials for the angular 
eigenvalues corresponding to unbound two-body systems 
(a s /b = —15.7 in fig. |l]). The external field in eq. (^) is 
negligible for small p and therefore the radial potential 
is negative when A + (37V - 4) (37V - 6)/4 < 0. Then 
genuinely bound states of negative energy are possible in 
our model without the confinement from the trap. 
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FIG. 2: Radial potentials Uo and Ui from eq. (|7() correspond- 
ing to the two lowest angular potentials for N = 20 and 
a s /b = -15.7 in fig. [j]. The trap length is bt/b = 1442. Also 
shown are the two lowest energies and the radial eigenfunc- 
tions fo and fo for the lowest radial potential. 

The radial potential for the lowest angular eigenvalue 
has a global minimum with U m - nl < at small p (« pb) 
separated by a barrier at intermediate p from a local 
minimum at p ~ p t = y/3N/2 b t . At large p the radial 
potential diverges due to the harmonic term, at small p 
due to the centrifugal term. The radial potential corre- 
sponding to the first excited angular state, also in fig. ||, 
coincides with the lowest radial potential at large p, but 
does not contain the attraction at small p (see fig. [I]) and 
diverges therefore to + infinity for small p without going 
through another minimum. 

The radial equation in eq. (^) has only one solution 
with negative energy E and the corresponding wave func- 
tion, shown in fig. ^, is located in the global minimum. 
The first of the infinitely many excited states in this po- 
tential is located in the local minimum at larger p created 
by the competition between the centrifugal barrier and 
the external harmonic oscillator potential. This excited 
state is usually referred to as the condensate, but the 
moderate attraction for the relatively small a s — —15.76 
already results in a lower lying many-body state about 
Pt/pb ~ 37 times smaller than the condensate. 

Large scattering length. At the threshold for binding 
of the two-body system the angular eigenvalue outside 
the interaction range must approach a negative constant 



Aoo, see fig. [j]. The scattering length is then in- 
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finitely large and the radial potential has the form 
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where £ 2 = -A^ - (3JV - 4)(37V - 6)/4 - 1/4 « lliV 2 /4 
for large N, where w — 57V 2 . Without the external p 2 
potential in eq. ( |lC)| ) the infinitely many radial solutions 
to eq. (p|) all behave like 



foo(p) = v^sin(|£|m(p/p sc )) 



(11) 



with some hyperradius scale p sc . The energies and sizes 
of the eigenstates labeled j are related by Q ^2| 
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With increasing quantum number these states become 
exponentially larger with exponentially smaller energies 
approaching zero. They originate from the constant A 
and the generic 1/p 2 potential in eq. Jlo|). They are 
many-body states, but not the embedded three-body Efi- 
mov cluster states which also arise from a 1/p 2 potential. 
However, both multi- and three-body Efimov states ap- 
pear for very large two-body scattering lengths. 

The external harmonic oscillator limits the possible 
number of these new states with E < 0. These negative 
energy states have to be located inside the trap and out- 
side the two-body potential. The number of nodes of 
allowed in this region equals the number of bound states 
Ne ~ |£|tt ln(p max /Pmin), where the first and last zero 
points then are at the end points of the interval. We 
obtain p max /Pmin ~ Pt/pb (see also fig. |) 



for b t < N\a s 



(13) 



where the quoted condition of validity is found from pt <C 
Pthj he. the plateau value must extend beyond the trap. 
The number of these available Efimov states therefore 
scales proportional to N (£ oc N) and logarithmically 
with the ratio of trap length and interaction range. 

We show in fig. || the radial potential for N = 20 and 
infinite scattering length corresponding to A,^ ~ —1340 
or £ 2 ~ 584, see fig. [j]. Using the estimate in eq. (|lj) 
we get Ne ~ 28 in good agreement with the computed 
30. The lower curve in the inset of fig. || is according 
to eq. ( |l2| ) a straight line of slope —1. The very lowest 
states deviate, since they "feel" non-constant A at small 
p, and the states close to E = deviate due to the in- 
fluence of the external potential. The energy spectrum 
becomes even denser above zero energy. For large posi- 
tive energies in the upper part of the inset the harmonic 
potential dominates and a straight line with slope +1 is 
obtained. The small positive energies are influenced by 
both external trap and interaction potentials. 
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Using eq. (@) we get 2(p 2 ) = {N - l)(rf 2 ) « 2(iV - 
l)(r 2 ). Even the most bound state, (p 2 ) 1 ^ 2 « 1366, has 
then a root mean square distance between two particles, 
(r^) 1 / 2 446, much larger than the interaction range. 
Also the root mean square radius (r 2 ) 1 ^ 2 ~ 316 is large. 
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FIG. 3: The lowest radial potential for N = 20, a s /b = oo, 
and bt/b = 1442. The horizontal lines indicate the 69 (30 
below zero) lowest energy eigenvalues and the inset shows 
their mean square hyperradii as a function of the absolute 
value of their energies relative to the values for the lowest 
state, i.e. 2mb 2 E /h 2 ~ -0.0147, \J(p 2 )o/b ~ 136. 

At the plateau, A ~ Aoo, £ 2 > and the radial poten- 
tial has no intermediate barrier. The centrifugal bar- 
rier can only be at p > p t \ n i.e. in the region where 
A ~ A n _ c , which provides the estimate of the barrier po- 
sition p ~ 2fcA 3 / 2 |a s |/3 > pth) see e q- @- A criterion of 
stability is then obtained by using A n _ c to estimate when 
the corresponding radial barrier disappears. This yields 
in analogy to Jl9| that the condensate will be (meta-) sta- 
ble if N\a s \/bt < 0.671. This critical value is from the 
Gross-Pitaevskii equation found to be about 0.6 fi p3| , 
and recently measured to 0.46 || for a 85 Rb-condensate. 

Stability is determined by decay through a barrier, like 
in fig. ||, when negative energy states and the condensate 
state in the second minimum are present. The tunnel- 
ing process then populates the states in the first mini- 
mum, i.e. many-body states of smaller extension and cor- 
respondingly larger density than the initial condensate. 
Competing with this macroscopic (all particles together) 
tunneling are the two- and three-particle recombination 
processes, which can occur from the initial condensate, 
or more likely from the collapsed many-body states at 
larger density and smaller p. Thus we have established 
a decay mechanism, i.e. macroscopic tunneling followed 
by recombination processes from the intermediate model 
independent states. 

An even more dramatic collapse could be initiated ex- 



perimentally by first creating a condensate correspond- 
ing to fig. ||, and then suddenly increasing the scattering 
length to that of fig. | §, 0, ||. The multi-particle 
Efimov states of fig. || would then quickly be populated, 
since the barrier is removed, and these states would sub- 
sequently leak out of the trap due to recombination. 

Conclusions. We use a new formulation designed to 
investigate two-body correlations in boson systems. The 
qualitative properties of Bose-Einstein condensates are 
explained by using the hyperspherical adiabatic expan- 
sion. Finite range and attractive potentials introduce 
new features compared to the mean-field approximation. 
Large scattering lengths are specifically treated. Many- 
body Efimov states appear for scattering lengths of the 
order of the trap length as model-independent states of 
much smaller extension than the condensate. Recom- 
bination through these intermediate states is suggested. 
This opens the possibility for realistic calculations of few- 
body recombination processes and for macroscopic col- 
lapse processes. 
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